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ABSTRACT 


A writ case assessment has been performed on the risk of an internal 
electrostatic discharge (ESO) on the Command Data System (CDS) multilayer orinted 
circuit board and an output power transformer module in. the oower subsystem 
cf the Galileo Orbitor. 

An estimate of the Jovian environment during the 35 hour orbit insertion 
was supplied by JPL and used as an input to calculate the electron transport 
into the Galileo components. A radiation shielding analysis computer code, 
CHARGE, calculated the electron transport deposition trapped in the 
anticipated sensitive areas of the multilayer board and transformer module. 

Based on these trapped charge calculations electric fields were calculated 
between the identified isolated areas and the spacecraft ground. These 
fields were then combined with the characteristic, dielectric properties 
of conductivity and dielectric strength of the materials to assess the 
risk of an ESD. In most cases, the deposition rate during the 35 hour 
encounter was neglected and was only included in the analysis 
of the results where the calculated fields approached the materials 

dielectric strength. 

The results of the assessment of ESD in the CDS multilayer printed circuit 
board indicate that the probabiMty of ESD in the FR4 is low. The probability 
of ESD in the components attached to the multilayer board, however, is 
uncertain based on a lack of prior experihental data. No analysis of any 
components were performed during this task. 


An experiment 1 : reconinended using one or two of the CDS prototype 
multilayer printed circuit boards to assess the risk of ESD in or on the 
multilayer board components. 

The results of the risk assessment of an ESD In the CDS transformer 
module Indicate that If construction of the transformer results In a 
nominal spacing between the core and the metal posts In the module the 
probability of ESD after a single orbit Is low. but the probability 
Increases as the spacing gets close to 5 mil. This risk can be reduced 
to zero by the addition of a grounding strap from the core to the space- 


craft ground. 
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INTRODUCTION 


High energy electrons and protons trapped within Jupiters magnetic fields 
and incident on the Galileo orbiter during its orbit insertion have the 
potential of causing interference and possible damage to the spacecraft. 

The work reported here is a worst case assessment of the risk of an 
electrostatic discharge (ESO) within the Galileo orbiter due to the 
accumilation of these energetic particles. 

Although the low energy «40KeV) electron charging process that produces 
external discharoes has received considerable attention over the past 
five to six years and a significant amount of analytic and experimental 
data has been generated, a comparable level of understanding of the ESD 
hazards resulting from the high energy charged particle spectrum has not 

been achieved. 

It was therefore, decided early in the program that a worst case assessment 
of a limited number of components of special concern would provide the best 
utilization of resources. This approach was used to most efficiently 
identify the magnitude of the problem and the need if any, for additional 

detail. The areas of needed detail were identified as either geometry 

definition or empirical or material data. 

The areas of special concern were identified at the kick off meeting at 
JPL. The areas were defined acconling to the current designs, as far as 
geometry and materials, with support from JPL and GE personnel in the 
orbiter program. This definition was then used in computer transport 
calculations at GE to determine the electric charge trapped during the orbit 
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Insertion. The worst case calculations assumed that the performance would 
be Independent of rate or Induced conductivity. The experience at 6E and 
Jaycor was then applied to assess the risk of the resulting generated 
fields leading to Internal ESD. 


II. TECHNICAL DISCUSSION 
2.1 Technical Approach 


2.1.1 Task Description 

The basic approach of the Internal ESD risk assessment was outlined at the 
kickoff meeting with JPL, 6E and Jaycor at JPL. The assessment was defined 
In seven tasks to track the progress through the effort. These tasks are 
described below and their Interrelation Is shown In Figure 1. 


Task 


Identify a limited number (2 - 3) of mission critical or 
sensitive areas of the spacecraft which In the opinion of 
the JPL and GE/Jaycor could serve as a first 
case assessment of the potential risks to the Galileo 
orbiter. 


This task will define most probable electron trajectory 
and the dielectric materials of greatest concern for 
potential Internal ESD to the Galileo Orbiter. 


Task 2 Confirm or update the Galileo radiation environment In 
both eneroy and temporal dependence to be used as the 
source of' Internal charging high energy electrons. 


Task 3 Establish a one dimensional geometric shielding and 

description for each of the scenarios Identified in Task i. 
This Usk will define all the geometries and materials 
parameters which affect the high energy electrons penetration 
into the dielectrics defined In Task 1. 


Task 4 Results of Task 2 and 3 will be Input Into the shielding 
code "Charge" to determine the resultant charge deposition 
as a function of depth and magnitude. 

Task 5 The output of Task 4 will be used to determine the electric 
f1^ degenerated Inside the dielectrics and will be compared 
to the existing materials data and ESD data base. 

Task 6 The computer calculations and material properties <^«ul ts of 

Task 5 will be compared to provide a risk assessment of Internal 
ESD or Identify whether additional materials property data 
Is needed or more refined computer calculations. 
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In the power subsystem resulted In a reassessment of the priority Ust. 

The following list describes the final priority list for Internal ESD case 
studies. 

Cise 1 The first three areas of Initial concern were combined 
Into a single case study of the electron charge 
deposited In the CDS multilayer board located In Bay A 
of the despun section where the electron trajectory of 
most Interest would be normal to the print circuit 
board. 

Case Z Consider the electron charge deposited In the core of a 
transformer In the power circuit In Bay B of the despun 
secti on . 

Case 3 Consider the electron charge deposited In the cable 
Insulation In the harness bundle between the spun and 
despun sections. 

2.1.3 Environment 

From data supplied by Nell Divine of JPL the electron {and proton) flux 
environment was obtained for L 4 16 (L < 12) along the proposed Galileo 
orbit Insertion. The flux was for electron and proton energies between 
63KeV and lOOMeV in logarithmicly spaced energy intervals where A log E 
■0.2. Figures 2 and 3 show the temporal dependence of the electron flux 
(electrons/cm^ sec) and fluence (electrons/cm^) for three particular 
energy levels. As a worse case initial assessment of the hazard from the 
deposited electron charge, the rate dependence was neglected. The total 
fluence as a function of energy was used from this Task 2 effort as an 
inout to the charge deposition computer code described in the nexv section. 
As part of this code input the total incident charge deposited during the 
nearly 35 hour encounter was assifned to be deposited in one hour. Table 
1 summarizes the total electron fluence (electrons/cm ) and assumed electron 
flux and differential electron flux used as input to the computer code 
CHARGE. Comparison of the assumed electron flux shown in Table 1 assuming 
all deposition occurs in one hour and the peak electron flux for the 
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TABLE 1 

ELEaRON ENVIRONMENT USED IN CHARGE CODE 


Energy 

(^in) 

(MeV) 

Fluence (N.) 

“e > E 
^ ^ ^MIN 

FLUX (J) 
(cm"^ sec'b 

DIFFERENTIAL ENERGY 
(dJ/dE) 

(on*^ sec'^MeV'^) 

100 

1.479 

E9 

4.108 

E5 

•6.102 

E4 

63 

8.128 

E9 

2.258 

E6 

-1.973 

E5 

40 

4.265 

E10 

1.184 

E7 

-1 .458 

E6 

25 

1.995 

Ell 

5.542 

E7 

-6.16S 

E6 

16 

5.754 

Ell 

1.598 

E8 

-1.708 

E7 

10 

1.122 

El 2 

3.117 

E8 

-3.215 

E/ 

6.3 

1.698 

El 2 

4.717 

E8 

-5.170 

E7 

4.0 

2.239 

El 2 

6.219 

E8 

-9.411 

E7 

2.5 

2.884 

El 2 

8.011 

E8 

-1 .808 

E8 

1.6 

3.802 

El 2 

1.056 

E9 

-4.173 

E8 

1.0 

5.138 

E12 

1.427 

E9 

-9.381 

E8 

.63 

7.079 

El 2 

1.966 

E9 

-2.250 

E9 

.4 

1.000 

El 3 

2.777 

E9 

-5.147 

E9 

.25 

1.412 

El 3 

3.922 

E9 

-1.260 

E10 

.16 

2.089 

El 3 

5.802 

E9 

-3.240 

ElO 

.10 

3.162 

E13 

8.783 

E9 

-8.730 

ElO 

.063 

5.138 

El 3 

1.427 

E10 

-1.371 

Ell 
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three energy values plotted In Figure 2 shows that the assumed value is 
higher than the peak flux by less than a factor of 10. 

Figure 4 shows the assumed electron flux given In Table 1 and the 
differential flux (electrons/cm‘ sec with energies between and 
calculated from Table 1. The differential flux Is the difference between 
fluxes at adjacent energy values. Also shown Is the agreement between 
the calculated ambient differential flux between energy levels and that 
calculated by the computer based upon the Input of the differential energy 
flux spectrum (dJ/dE) given In Table 1. 

2.1.4 Electron Transport Code 

The results of Tasks 2 and 3 were used as inputs to the radiation transport 
code CHARGE. The original intent of CHARGE was the Space Radiation analysis. 

The GE Space Division code package is a modification of one developed at 
Douglas Aircraft Co. in the 1960's.^ This code package. Is reasonably 
economical in terms of input and computer expense. As one of its outputs 
the code supplies the incident and exit current densities as a function of 
energy interval in addition to the total fluence. The differences between 
these values were used to calculate the charge deposited in a particular 
layer. The following discussion describes the operation of CHARGE. 

The CHARGE code computes response rates and flux spectra behind a multi- 
layered spherical or planar shield exposed to isotropic fluxes of electrons, 
protons, and heavy charged particles. The doses or other responses to 
electron, primary proton, heavy particle, electron Bremsstrahlung secondary 
proton, and secondary neutron radiations are calculated as a function of 
penetration Into the shield; the materials of each layer may be mixtures of 
elements contained In the library or supplied by the user. The primary 
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results consist of flux responses as a function of penetration. 

The ambient electron, proton, and heavy particle spectra may be specified 
as differential spectra In tabular or functional form. These Incident 
charged particle spectra, which may be viewed as essentially continuous 
functions of particle kinetic energy, are divided Into energy bands or 
groups, the number and spacing of which are controlled by Input data. Each 
of these groups may be completely described by Us upper and lower energy 
limits and the energy variation of the spectrum within the group. 

The variation of the group boundary energies and group spectra as a function 
of shield penetration uniquely determines charged particle dose rates and 
secondary particle production rates. Therefore, charged particle shielding 
calculations reduce to the Integration of the equations which express the 
variation of particle energy with distance In the shield material ("range- 
energy" equation). 

As charged particles penetrate matter, they Interact primarily with the 
orbital electrons of the medium. Incident electrons can undergo quite large 
deflections in their billiard ball -like collisions with the orbital electrons. 
In addition to producing electromagnetic radiation (bremsstrahlung X-rays) by 
the deceleration of these charges, the Individual large angle collisions 
can cause large energy losses and deflections from the original electron 
path. In order to account for these competing processes, CHARGE uses the 
j basic range-energy relation, normally modified by applying electron 

^ transmission factors derived from Monte Carlo calculations. The capability 

of not using transmission factors Is contained, however. 

Because of the limited number of materials defined in the computer library 
and the relative similarity In the elemental breakdown of mylar. Teflon and 
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soli thane, these materials were grouped together under a cownon description 
of polyethylene which existed In the code library. 

2.2 CDS Multilayer PC Board 
2.2.1 Material Definition 

The first case studied was the electron charging In the 10 layer multi- 
layer board Inside the COS In Bay A. The trajectory of the electron was 
defined to be normal to the CDS board. As such, the electron path was 
defined as passing through a multilayer Insulation blanket, the Bay A 
outboard cover, the board’s conformal coating and the multiple layers of 
the board's dielectric. 

The description of the multilayer Insulation blanket shown In Figure 5 
was supplied by JPL. The description of the outboard cover on Bay A 
parallel to the PC board was taken from Drwg. No. 10091629, Sheet 2. The 
construction detail of the multilayer board was taken from JPL Drwg. No. 
10064876, Sheet 1 and discussions with JPL during the kickoff meeting. The 
boards were described Lo be 10 to 16 layers (typically 10 layers) of 
materials with a 20 mil conformal coating of Soli thane 113-300. The 
dielectric spacing Is 8 mil between conductor layers. The conductor strip 
line Is 1.4 mil thick and 8 mil wide. The typical lay up Is shown In 
Figure 6. Based upon these descriptions the material profile through which 
the electron trajectory passes Is given in Table 2. 


TABLE 2 


LAYER NO. 

DESCRIPTION 

THICKNESS 

IN. (mm) 

1 

Black Paint 

0.0015 

(0.038) 

2 

Tef 1 on 

0.002 

(0.050) 


Mylar 

0.001 

(0.025) 


15 Layers Mylar (1/4 mil) 

0.00375 

0.095) 


Mylar 

0.001 

(0.025) 

3 

Tantalum 

0.006 

(0.152) 

4 

A1 umlnum 

0.020 

(0.508) 
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.001 HYLAR, ALUMINIZED 


INNER BLANKET 
OUTER BLANKET 



15 LAYERS .00025 MYLAR 
ALUMINIZED BOTH SIDES 

16 LAYERS .007 DACROM NET 

.001 MYLAR, ALUMINIZED 
BOTH SIDES 
.002 TEFLON ALUMINUM 
BOTH SIDES 

.0015 CONDUCTIVE BLACK 


TYPICAL RPM THERMAL 
BLANKET CROSS SECTION 


SCALE: NONE 


FIGURE 5 THERMAL BLANKET 



INSULATOR 


FIGURE 6 MULTILAYER BOARD 
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LAYER 


DESCRIPTION 


THICKNESS IN. (m) 

.5 Soli thane 0.020 (0.508) 

5 to 15 0.008 (0.203) 

The vapor deposited aluminum (VOA) on all the teflon and mylar films was neglected 
In these calculations because of their negligible cross-sectlor*"* area for 
electron capture. The thickness of the tantalum shield was chosen from Orwg. 

No. 10091629. It Is the minimum shield thickness defined over the Bay A cover. 
The alumlnun thickness was assuned to be 0.06 Inches but from later discussions 
with JPL this value may be low by 0.02 Inches. Per our discussion with JPL the 
so11t:;ane thickness was assumed to 0.02 Inches. The next ten layers were assuned 
to be FRA glass epoxy each 0.008 Inches thick. The Inclusion of a discontinous 
copper layer about 0.0014 Inches thick between layers of FRA was not Included 
In the analysis because of the lack of a good definition for any "typical 
profile of conductor through the board for a given cross section. This analysis 
allowed a worst case construction assuming a cross section which contained no 
Intermediate ground planes for charge leakage. 

The FRA was described by an elemental breakdown of the largest constituents 
using Vendor data according to the following table: 

TABLE 3 

ASSUMED ELEMENTAL COMPOSITION OF FRA 


Oxygen 

33% 

Carbon (+ Nitrogen) 

31% 

Aluminum (+ Silicon) 

19% 

Krypton (Bromine) 

9% 

Calcium 

8 % 


lA 


The elements In {>arentheses were also In the material but were described 
or combined with the first element because of lack of radiation shield data 
un these materials. The density of FR4 was assumed to be 1.93g/cc. 

2.2.2 Charge Deposition Profile 

Since the charge deposition was modeled In a planar geometry the code used 
one half the Isotropic flux given by the dashed line In Figure 4. For each 
layer described In Table 2 the CHARGE code calculated an Incident and exit 
total fluence and an electron flux as a function of renormalized energy 
Intervals. Figure 7 shows the Incident and exit flux as a function of the 
electron energy for the first FR4 layer. Table 4 sums the total electron 
fluence Incident and penetrating each of the material layers. Differences 
In these values represent the trapped charge within each metal or dielectric 
layer assuming no leakage due to Its electrical conductivity. 

With the exception of the first FR4 layer the charge density is relatively 

10 2 

uniform with about 2.6 X 10 electrons/cm per 8 mils or .2mm of multilayer 

board. This gives a charge density of 2 X 10”^ Coul/cm^. The total trapped 

charge In the 100 mil thick composite board and sollthane layer Is nearly 3 X 10^^ 

2 

electrons/cm . 

2.2.3 Internal Field Calculations 

The results of the environment transport calculations In terms of the trapped 
charge density were used by Jaycor as Inputs to the PRECH6 computer code (See 
Appendix) to determine the electric fields within the multilayer PC board. 

Since the electric fields are a function of the relative position and density 
of the conductors In the multilayer PC board which vary throughout the actual 
PC board. It was not practical to analyze all configurations that could exist 
on an operational board. Therefore, two worse case geometrical configurations 
were assumed, one with a single grounded conductor on the surface of the 10 layered 
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TABLE ^ 




INCIDENT ELECTRON DENSITY 

DEPOSITED CHARGE 

LAYER 

DESCRIPTION 

{e/cm2l 

(e/cm2) 

1 

Carbon 

2.569 

El 3 

5.857 

El 2 

2 

Blanket 

1.983 

El 3 

9.S43 

El 2 

3 

Tantalum 

9.889 

El 2 

7.860 

E12 

4 

A1 

2.029 

El 2 

6.815 

Ell 

5 

Sol 1 thane 

1.347 

El 2 

4.50 

ElO 

6 

FR4 

1.302 

El 2 

3.24 

ElO 

7 

II 

1.270 

El 2 

2.63 

ElO 

8 

II 

1.244 

El 2 

2.63 

ElO 

9 

II 

1.217 

El 2 

2.63 

ElO 

10 

M 

1.1 91 

El 2 

2.59 

ElO 

11 

n 

1.165 

El 2 

2.52 

ElO 

12 

II 

1.140 

El 2 

2.63 

ElO 

13 

II 

1.114 

El 2 

2.41 

ElO 

14 

11 

1.089 

El 2 

2.34 

ElO 

15 

II 

1.066 

E12 

2.23 

ElO 

16 


1.044 

El 2 

— 
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bo.rd o« «iM th« 9 «und.d conductor burlod In t(>. cnt.r of th. depth of 
the 10 l.yer hoerd. Since the electric fields neer the conductor strips 
increese es their separetlon betueen conductors Increeses, e uorst cese separation 
of one inch In the plane of the laultlleyer board «as assuaed. mile this Is 
«,t representative of the najorlty of the board, discharges In any portion of 
the board could couple Into other circuit elenents on the board. In addition 
the calculated electric fields can be fairly -ell scaled linearly es a function 
of conductor separation and thus the results can be applied to areas with 

closer spacings. 


The configuration uas modeled In a toroidal shape as shown In Figure 8. since 
PhECHG reoulres an axially sy^etrlc gemaetry. Figure 8. represents the surface 
conductor configuration while Figure 8b represents the burled conductor. The 
width of the dielectric was taken to be equal to one Inch, per the above 
discussion. The radius of the dielectric circle was made large compared to th. 
dielectric width so that th. electric fields calculated for these configurations 
sr. essentially Identical to those that would be calculated If the dielectric 
"tube” were unfolded Into a straight line. The E field on one side of the torol 
due to electric charge on the other side 1$ down by a factor of (r^/r^) a 0.01 
fr^ th. field due to the charge on th. same side. The conductor was assma.d 
to be 1.4 mils thick and 10 mils wide which Is th. anticipated dimensions for 
the flight boards. The spatial grid for which th. fields were computed was 
on th. order of 0.4 mil near th. conductor and became progressively larger at 

distances away from the conductor. 


As on. would expect, th. largest electric fields occurred near th. edges of 
th. conductor. Neglecting the conductivity of th. FR4. th. largest elKtrlc 
field was "Slculated to be 1.3 X loS/.11 (5 X lo‘v/cm). When th. conductor was 
burled at mid depth In the dielectric and about 4.3 X 10 V/mll (1.7 X 10 V/cm) 
when th. conductor was on th. surface. If the corners of th. metal strip 
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Top view 
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FRA (1 In. wide) 

Metal Strip (10 nil wide) 


Side View 
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1.6 mil 

80 mil • i 


NOTE: This drawing Is not to 
scale 
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FIGURE 8 TWO CONFIGURATIONS FOR CALCULATION OF E FIELD INSIDE FRA BOARD 
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were not perfectly sham as assumed In the code, these two fields would be 
closer to about 8,9 X lO^V/mll (3.5 X 10®V/cm) and 3 X lo\/mil (1.2 X lo\/cm) 
respectively, which are approximately the average fields over the surface of 
the conductor. 

Hhen the conductivity of the FR4 Is Included, a calculation of the conducted 

current In the plane of the FRA can be performed and compared to the Incident 

r 

current density. Using an average electric field of 2 X 10''V/cm, which was 

• 1 3 

calculated assuming no conductivity, and a value of 2 X 10' (^cm)' for the 

conductivity, one finds that the lateral conducted current Is 
J(A/cm^) ■ E(V/m) ♦(!( (A-cm)'^) 

• 2 X 10® • 2 X 10’^^ 

■ 4 X 10"®A/cm^ 

To produce this current In the plane of the board, the Incident flux normal 

• 9 2 

to the board would have to be 9 X 10“ A/cm . This value Is roughly 700 times 

2 

greater than the maximum anticipated Incident flux of 13pA/cm . Thus with 

a value for the conductivity of 2 X 10“^^ (o.-cm)"^ the maximum fields Inside 

the multilayer board will In reality be roughly 700 times smaller than those 

calculated with zero conductivity. The maximum fields will then be on the 

order of 13V/m11. This Is close to the value calculated assuming a close 

packed density of conductors between the layers of FRA and using the calculated 

*•8 2 

charge density of 2.6 X 10 electrons/cm In the 8 mil FRA layer. 

2.2.4 Risk Assessment 

An assessment of the risk of an ESO occurring In the COS (Bay A) mulltlayer 
PC board can be performed In two ways. 1) The results of the environment 
transport calculations can be used In analytical computer codes to calculate 
electrical stresses In the multilayer board, and compared to the material 
electrical stress characterizations. 2) The results of the environment 

transport calculations In terms of trapped charge In the multilayer PC 
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board can bt coaipared to tx1 sting axparlmental data on discharges in printed 
circuit boards. Both approaches are treated here. 

Coeiparison of the calculated rate of charge trapping { the flux of trapped 
electrons) in the Jovian environment to Figure 9 is one means of assessinr 
the risk of ESO in the multilayer board. From Table 4, roughly 3 X 10 
electrons/cm^ will be deposited during a 35 hour Jupiter encounter. The 
deposition time for the calculation was assumed to be one hour, which implies 
a deposition rate of 8.3 X 10^ electrons/cm^ sec (13pA/cm^). This flux, however, 
is artifically high since the 3 X 10^^ electrons/cm^ will be deposited 
continually during the 35 hour encounter with the Jovian electron environment. 

The actual flux will be low at initiation of the orbit insertion, will peak 
near closest encounter and will again decrease as the Galileo moves away from 
Jupiter. Thus, the majority of the trapped charge will be accimulated during 
a time between the one hour assumed for calculational purposes and the actual 
35 hour encounter time. The value of 8.3 X 10^ electrons/cm sec is certainly 
an upper limit on the flux. A lower limit for assessment of the risk of ESD 
in the multilayer board should be roughly one-tenth this value or 8.3 X 10 
electrons/cm^ sec (1.3 pA/cm^). Comparison of these values to Figure 9 indicates 
that discharges could be expected to occur roughly once every 850 sec (14 minutes) 
at 8.3 X 10^ electrons/cm^ sec (13pA/cm^) and once every 17 hours (extra- 
polated from inverse square fit to the data) at 8.3 X 10® electrons/cm sec 
(1.3pA/cm*). Between these two limits based on a comparison of Table 1 and 
Figure 2 it is anticipated that several discharges could occur during an hour. 

One the other hand, based upon electron transport calculations and analytical 
computer code predictions of the electrical stress in the FR4 material the 
electric fields predicted in the material are very small compared to the 
dielectric strength of FRA of roughly 500 volts/mil. 
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At this tine only United data exists on the charging of Internal spacecraft 
conponents by high energy, penetrating electrons. The only existing data on 
the charging of printed circuit boards was obtained by Jaycor during experinents 
perfomed for the Air Force Weapons Laboratory (Ref. 2, 3, 4 ). The energy of 
the electrons used In these PC boards tests was about 2S0KeV and resulted In 
a charge which was distributed throughout the first 33 mils of the PC board. 
Therefore, these results hold a large significance to the risk assessnent. 


Briefly stated, the results of these experinents were: 

(T) ESD can be Induced In or on satellite PC boards by high 
energy electrons. 

(2) During ESD In or on PC boards currents and voltages are 
often Injected into electronic components which are' well 
above those necessary to cause upsot In nost types of 
electronic components and sufficiently large In some cases 
to cause burnout. 

(3) Both the magnitude and frequency of the oischarges de> 
crease as the flux of Incident electrons (rate of charge 
trapping) decreases. Figure 9 shows the time between 
discharges (the Inverse of the discharge frequency) as a 
function of Incident electron flux. Note that as the 
electron flux decreases, the time between discharges 1n> 
creases roughly quadratically. 

(4) Integrating the Incident flux over the time from the 
start of the electron exposure to the time of the first 
ESD gives a threshold electron charge layer density for 
ESD's of 1.5 X loll electons/cm^. 

Although the charge calculated for the CDS multilayer PC board In the Jovian 
environment will be distributed throughout the entire 90 mils of the 
multilayer board, the resultant electric fields will not differ by more than 
10 to 20t from the fields present In the experiments described above where 
the charge was distributed through only the first 33 mils (<^1/2 the thickness) 
of the PC board. Thus, a comparison of the calculated charge density In the 
COS board to the ESD threshold charge density Implied from the experimental 
data Is valid. Comparison of the experimentally determined ESD threshold charge 
layer density of 1.5 X 10^^ rlactrons/cm^ to the roughly 3 X 10^^ electrons/cm^ 


22 


In Table 4, which will be trapped In the so11thane/10 layer FR4 PC sandwich 
Is Indicative of a high probability of ESD's occurring In or on the CDS multi* 
layer PC board. 


In view of the electric field calculations when consideration Is given to 

the relatively high conductivity of FR4, the probability of £SD In the CDS multi* 

layer board Itself Is very small. The apparent discrepance between the 

experimental results which Indicate that ESD In the test PC board will occur 
2 

at rates of 13pA/cm and the electric field calculations results for the 

multilayer FR4 board could be due to one of two things. 

The printed circuit board for which ESD data does exist had a 
conductivity much lower than that of FR4. The test printed 
circuit board was fiberglass. Standard varltles of printed 
circuit board fiberglass have conductivities that range from 
10-13 to 10- “*6 (.n.*cm)-^ 

The ESD observed during the printed circuit board tests did not 
occur In the printed circuit board Itself, but In the components 
on the board. In particular, one could Imagine that the ESD 
occurred In capacitors on the board. This Issue has to date not 
been explored. 


The probability of an ESD in the components (possibly capacitors), therefore, 
should be addressed before drawing a conclusion regarding the probability 
of an ESD In the CDS electronics. 


2.3 Power Transformer Module 
2.3.1 Material Definition 

The second case studied was the electron transport through the Bay B. 


outboard cover and thermal control blanket Into a power transformer module. 
The concern for Internal ESD In the transformer module was associated with 
the electrically floating core. The anticipated thermal control cover of 
Bay B Is a combination of thermal louvers and multilayer Insulation (HLI) 
blanket. The structural support In the Bay B cover required to support the 
louver and provide other particulate protection was specified as a 192 mil 
thick wall of aluminum. Therefore, the worst case for electron charging 
Inside Bay B would be to consider an electron trajectory through the MLI 
blanket similar to the profile used for the multilayer board. For this case 
the thickness of the aluminum cover was adjusted to 80 mils to reflect the 
latest Information on the shielding description. 

As In the previous case the worst case orientation of the transformer module 
was assumed to be located under the area of the Bay B cover with the minimum 
of tantalum shielding. 

Initial drawings from the Voyager program of the transformer Indicated 
Insulation tape between windings. Further discussions with JPL and Jeff 
Benham of GE about the geometric details of the transformer core and windings 
resulted In the analyses being redirected to consider the anticipated new 
design for the Galileo output power transformer. Transformer dimensions 
and materials were obtained from the preliminary transformer Drwo. No. 
(GE-47E248720). Figure 10 shows the orientation of transformer core in 
the encapsulate module. 

The electron transport through the MLI blanket and metal cover Is similar 
to Table 2 with layer 4 becoming 80 mils thick. The normal electron path 
Is then through a nominal 2" mil coating of soli thane over the module. 







through *31 mil fiberglass mounting board and into the potting compound 
of soli thane or SHRO above the transformer. To calculate the amount of 
charge trapped in the core, two cases were treated for the electron path to 
the core. 1) The electrons traversed parallel to the edge of the windings 
through 287 mils of potting compound between the fiberglass connector plate 
and the area of the core marked as A in Figure 10. 2) The electrons 
traversed through 107 nils of potting compound between the fiberglass and 
the windings, through 11 mils of fibrous tape, through 6 layers of magnetic 
wire, through 15 mils of paper press board bobbin and into the core. The 
ratio of these two areas is about 3.08cm‘ (area A) to 2.90cm (area B). These 
two profiles are described in detail in Table 5. 


TABLE 5 


layer no. DESCRIPTION THICKNESS IN (mm) 


1 

Black Paint 

0.0015 

(0.038) 

2 

ML I 

0.0077 

(0.195) 

3 

Ta 

0.006 

(0.152) 

4 

A1 

0.080 

(2.032) 

5 

Sol i thane 

0.020 

(0.508) 

6 

FRA 

0.031 

(0.787) 

7 

Soli thane 

0.107 

(2.718) 

Case 1 

8 

Fibermat Tape 

0.011 

(0.279) 

9 

6 Layers of Wire 

0.154 

(3.912) 

10 

Paper Pressboard 

0.015 

(0.380) 

11 

Core 

1.000 

(25,400) 

Case 2 

8 

Soli thane 

0.180 

(4.570) 

9 

Core 

1.000 

(25.400) 
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2.3.2 Charge Deposition Profile 

Since the charge deposition Into the power transformer was modelled In a 
planar geometry the code used one half of the Isotropic flux given In Figure 
4 for the Incident flux on Bay B. For each layer described In Table 5 the 
CHARGE code calculated an incident and penetrating total fluence and 
electron flux as a function of a set of renormalized energy Intervals. 

Since the current design of the output power transformer In the power 
subsystem 1o Bay B of the Galileo orbitor has the core effectively Isolated 
from the spacecraft ground, the analysis concentrated on the charae deposited 
and trapped in the transformer core. Table 6 sums the total electron 
fluence Incident and penetrating each of the material layers. Differences 
In these values represent the trapped charge within each metal or dialectic 
layer assuming no leakage due to electrical conductivity. 

2.3.3 Internal Field Calculations 

The potential, V, which will develop across the transformer core and the 
spacecraft is given by V • J where Q is the amount of charge trapped in the 
conductive core and C is the capacitance between the core and the nearby 
metal. From Figure 10 it can be seen that the closest metal to the trans- 
former core is the transformer windings and the metal posts which are used 
to package the transformer. Since the widnings are much closer to the core 
than the four metal posts and since the surface area projected to the core 
by the windings is much larger than that projected by the posts, the 
capacitance of the core to the spacecraft will be dominated by the presence 
of the windings. This capacitance was approximated by considering each 
portion of the core covered by windings as a parallel plate capacitor. The 
capacitance Is given by C • ^ where d Is the spacing between the winding 
and the core, A is the surface area of the winding and e is the dielectric 
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TABLE 6 


INCIDENT ELECTRON DEPOSITED CHARGE 

LAYER DESCRIPTION DENSITY (e/ctn^) (e/cm^) 




Case 1 

Case 

2 

Case 1 

Case 2 

1 

Carbon Paint 


2.569E 13 



S.857E 12 

2 

Blanket 


1.983E 13 



9.941E 12 

3 

Tantalum 


9.889E 12 



7.860E 12 

4 

Aluminum 


2.029E 12 



8.15 E 11 

5 

Sol 1 thane 


1.214E 12 



3.80 E 10 

6 

FRA 


1.176E 12 



1.07 E 11 

7 

Potting 


1.069E 12 


1.68711 

1.557E 

8 

Winding 

9.003E 

11 


1.375E 

11 


Potting 


9.1 33E 

11 



9 

Insulation 

7.628E 

11 


2.100E 

9 

10 

Winding 

7.607E 

11 


5.540E 

11 

11 

Press board 

2.067E 

11 


2.500E 

9 

12 

Core 

2.042E 

11 7.452E 

11 

2.039E 

11 7.440E 

13 


2.803E 

8 1.194E 

9 


• • 
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constant of the paper pressboard and wire Insulation between the core and 
the winding. Using the dimens if^ns shown in Figure 10, the capacitance was 
calculated to be 2.6 x 10“^*^f, assuming a value of 4 for the dielectric 
constant of the pressboard/wire insulation combination. The relative 
dielectric constant for most cellulose based compounds ranges from 2 to 
6 at IMHz. From the results of the environment transport calculations in 
Table 6, 4.6 x 10’^ Coul are anticipated to be trapped in the core during one 
35 hour encounter with the Jovian electron environment. Ignoring the con- 
ductivity, the potential which will develop between the core and the windings 

will be 

4.6 X 10*^ coul 

2.6 X 10"’° f 

« 1,769V 

The electric field across the pressboard/wire insulation (total thickness 

17 miles) will be 
1,769 


• 104V/mil 

Subsequent encounters with the Jovian electron environment will possibly 
increase the electric field if no conductivity is assumed for the dielectric 
between the core and winding. However, at this point, one must consider 
conductivity effects since the time between consecutive orbits will be greater 
than 35 hours. If a conductivity of 10"^® (Acm)"^ is assumed, the leakage 
current after a single orbit using a maximum field of lOOV/mil, can be 

calculated to be 

I - E- • A 
» 4 X 10‘^°A 

Since this value is roughly 3 times larger than the anticipated peak current 
stopped in the core of 1.3 X 10’^° amperes, the peak fields will be reduced 
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by the conductivity by ot loast a factor of three, i.e. the fields resulting 
from a single encounter will be on the order of 30V/mil. The effects on the 
electric field of multiple orbits can be assessed by calculating a dielectric 
relaxation time which is give by 

'y m e 

• 308 sec 

Since this relaxation time is very short compared to the times between encounters 
with the intense electron fluxes (>35 hours) the effect of multiple encounters 
will be negligible. 

Recognize however, that if the conductivity were on the order of 10 ^®(.Acm) ^ 
then the peak fields would be on the order of that calculated assuming no 
conductivity. This point is significant because the value of 10 ^^(.r».an) ^ 
is representative of this type of material, however the physical construction 
of the tranformer may result in the actual conductivity between the core and 
the windings being lower due to poor surface contact, etc. 

Additionally, the possibility of the flight geometry varying from that shown 
in Figure 10 must be considered. In particular, if the transformer core is 
not positioned synmetrically with respect to the four metal posts and is 
closer to one post than to the others, the capacitance of the core to the 
spacecraft will be altered as will the spacing across which the core potential 
will have to be stood off. From Figure 10 the nominal spacing between the 
core and the metal post was calculated as 37 mils. The capacitance of the 
post to the core can be calculated as the capacitance of a wire above a 
ground plane which is given by 

C • 

{f }j j 
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f • ; 1 

i ■ : ; 

where L Is the length of the wire, h 1s the height of the center of the 
) wire above the ground plane and a Is the radius of the wire. If It Is assumed 

that the core Is packaged such that one post Is only 5 mils from tne core, 

' then the capacitance of the core to the post can be calculated to be 

K 1.6 X 10'^’ f 

assunlng a relative dielectric constant of 3.5. Since this Is much smaller 
than the capacitance of the core to the windings, the capacitance of the 
^ core to the spacecraft Is still dominated by the capacitance of the core to 

I the windings. The potential which develops between the core and the space- 

craft will have to be stood off by only 5 mils of the potting epoxy which 
corresponds to an electric field of 340V/m11. Consideration of the effect 
of the conductivity of the pressboard could possibly reduce this field by 
a factor of 3 as discussed above, thus the electric field may be only 
‘ 113V/m11; however the same consideration discussed above regarding the 

accuracy with which the conductivity of the pressboard Is known must be kept 
In mind. 

2.3.4 Risk Assessment 

The dielectric strength of most acetate-based materials ranges from 200 to 
600V/m11. Thus It appears that the risk of ESD due to charge trapped in a 
transformer core on one encounter with the Jovian electron envlroranent Is 
relatively small, since the field across the pressboard will be only about 
30 to lOOV/mil. 

In the case of an assynmetrical plac*r.£nt of the core In the transformer 
module where the core Is close to one of the metal posts as discussed In 
the previous section, the potential which develops between the core and the 
spacecraft will now have to be stood off by 5 mils of the potting epoxy which 
corresponds to an electric field of 113 to 340V/m11 depending upon the 
conductivity value assumed. The dielectric strength of common epoxies Is 
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betMcen 350 to SOOV/mll. Thus, tf the fabrication control specifications 
are such that the core can conceivably be within 5 mils of the mounting 
post, there is a risk of ESO's occurring from the transformer core to the 
metal posts during a single encounter. Although such an arc will not 
Involve direct current injection into the transformer windings, the arc will 
couple a current into the windings. If the 1769 V potential is reduced 
to zero, i. }. lOOX of the charge trapped in the core arcs to the metal post, 
the current induced in the windings can be calculated by 

I - Cy ^ 

" Ht 

If it Is assumed that the arc has a risetime of lOnsec 
I » 2.6 X 10‘^° • 1769 

■ 46 A 

In suninary, if the trans'former is fabricated as indicated in Figure 10, the 
probability of an ESD during a single encounter with the Jovian electron 
enviroronent is relatively low. If the configuration control on fabrication 
of the transformers is such that the core may be positioned within several 
mils of the metal posts, the probability of an ESD is significantly increased. 
If a path for draining the charge trapped in the transformer core to the 
spacecraft ground were provided, the risk of ESD from the transformer core 
would be effectively reduced to zero. 

3.0 CONCLUSIONS AND RECOftiENDATIONS 

A worst case risk assessment has been performed on the hazard of an internal 
ESD on a COS multilayer printed circuit board and an output power transformer 
core in the power subsystem of the Galileo Orbiter. 


32 


s 


An estfnate of the Jovian environment during the 35 hour orbit Insertion 
was supplied by JPL and used as an Input to calculate the electron transport 
into the Galileo components. A radiation shielding analysis computer code, 
CHARGE calculated the electron transport deposition trapped In the 
anticipated sensitive areas of the multilayer board and transformer module. 

Based on these trapped charge calculations electric fields were calculated 
between the Identified Isolated areas and the spacecraft ground. These 
fields were then combined with the characteristic dielectric properties 
of conductivity and dielectric strength of the materials to assess the risk 
of ESO. In most cases, the deposition rate during the 35 hour encounter 
was neglected and was qualitatively Included In the analysis of the results 
where the calculated fields approached the materials dielectric strength. 

3.1 Multilayer Board 

The results of the risk assessment of ESD in the CDS multilayer printed 
circuit board indicate that 

1. The probability of ESO in the FRA is low 

2. The probability of ESD in the components attached to the 
multilayer board is uncertain. 

Since the existing experimental data on printed circuit board ESD does not 
distinguish between ESD in the printed circuit board or in the components 
on the board and since the data does indicate that potentially damaging ESO 
did occur either in or on the printed circuit board, the risk of ESD in the 
Galileo CDS multilayer printed circuit board components should be examined 
further. Two methods of addressing this Issue are by analysis and by 
experiment. An analysis could make use of the previously performed environment 
transport calculations to calculate the electric fields inside selected 
components. An experiment would require exposure In a vacuum of a represent- 
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•tiv® multlliyer printed circuit bosrd to hlph energy electrons (200 - 
300KeV and possibly higher) at realistic fluxes, vVfth (heasurements of ESD 
currents and voltages at select locations. Both approaches have advantages 
and disadvantages. An analysis Is In general less expensive, but hOMever, 

****'^^^* generally limited to the availability of accurate materials 
property data. An experiment Is In general more costly, however, the results 
are generally more conclusive and convincing. 

The recownended approach Is to perform an experiment using one or two of 
the CDS prototype multilayer printed circuit boards. The experiment should 
be relatively straightforward and could possibly be performed using an 
electron source available at the Jet Propulsion Laboratory. 

3.2 Power Transformer Module 

The results of the risk assessment of ESD In the CDS transformer cores 
Indicate that 

1. If duality control on construction of the transformers results 
In a nominal spacing between the core and the metal posts on the 
order of 37 mil, the probability of ESD during a single encounter 
with Jovian electron environment Is low. 

2. If quality control on construction of the transformer Is such that 
the spacing between the transformer core and the metal posts may 
be as small as 5 mils, the probability of ESD during a single 
encounter Is significantly Increased. 

3. An ESD In the transformer couTd result In coupled currents In the 
Transformer on the order of 46A with pulse widths on the order of 
20ns. If the Injection of currents on the order of 46A for times 
on the order of 20ns Into the transformer windings are judged to 

be potentially deleterious to the functionino of the Galileo space- 
craft, the probability of ESD due to charged trapped In the 
transformer core can be reduced to aero by the addition of a grounding 
strap from the core to the spacecraft ground. Since the highest 
rate at which charge will be trapped In the core corresponds to a 
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egrrent 1t$» tMi. 130pA. the 9rogtKl1n9 ttrep eouU 

be e«trii»ely 1l9ht wi9ht. On the hesis of the resulte o t. 
H$k etsetsMht. the eddltlon of $»ch e 9round109 ttrep It 
rtcoiwended. If ESO induced currents in the traisfoneer wind ng 
on the order of *6A ere Judped poUntHUy hetardout to trtttlon 

performance. 

«. The pr1«ry uncerUintlet In thU ritk ettetsment retult fr» a 
lack of dielectric propertlet data and a firm knowledge of the 
transformer tolerances. 
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APPENDIX 


THE PRECBG CODE 

The PRECHG CODE vas developed io 1978 by JAYCOR as an analytical tool for 
the calculation of electric field distributions around cyllc :^rlcally aysnetrlc 
charged bodies consisting of conductors and/or dielectrics. The PRECHG code 
is essentially a PolBson*'Solver which computes the steady state electric 
fields on a finite-difference mesh grid Inside a cylindrical, conducting can 
due to a cylindrically-synnertlc distribution of charge on dielectric or 
conducting bodies inside the tan. The outer conducting can can represent 
either a real outer conductor or a boundary at Inifnlty if the can dimensions 
are made large compared to the internal bodies. The charge distribution 
can either be fixed in space, as in a dielectric, or it can reside in internal 
conducting bodies or in the conducting can walls. The charge in any 
conductor Is automatically adjusted by the code over the eurface of the 
conducting object to make the surface equlpotentlal. The total net charge 
inside the can Is always zero. In other words, if electrons are trapped in 
a dielectric or a conducting body, there is an Identical amount of positive 
Image charge somewhere else inside the can. 

A calculation with the PRECHG code is started by specifying an initial 
distribution of electric fields (E) along some, or all, of the mesh grid 
lines Inside the conducting can. This distribution of electric flux lines 
specifies the charge distribution of the problem. A positive and negative 
amount of charge ^Q) is assumed to be located at the two ends of the 
flux line. The magnitude of each 6.Q is the integral of the flux over the 
cross sectional area (Aa) that Is associated with the flux line. 
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^ • E* ^A, 

wlMte t Is the psrmlttlvity of the aedlum in which the flux line exists. 

If one end of the flux line starts at a conductor, A Q is pert of the total 
charge on that conductor. If the end of a flux line is not on a conductor, 
that AQ is assuned to be fixed at the position in space. 

This initial distribution of flux lines is somewhat arbitrary because there 
is an Infinite number of distributions that could be chosen to produce the 
desired charge distribution. In principle, any flux distribution that is 
consistent with the desired charge distribution could be used but usually 
the simplest arrangement is chosen. This distribution of flux lines satisfies 
conservation of charge (and electric flux) but it normally wiU not satlfy 
*7 X E • C, which is necessary condition for satisfying Poisson's equation. 
The primary computational process in the ?RECHG code is an iteration of the 
magnitudes of V X E every mesh point Inside the conducting can in such a 
way that 7 X E tends toward zero everywhere while conserving charge at 
each mesh station and in each conducting body. The iteration process is 
terminated when the largest magnitude of V X E at any place in the volume is 
below some preselected value or when the values of E are no longer changing 
significantly with successive iterations. This final distribution of E is 
the solution of Poisson's equation for the chosen distribution of charge, 
within the accuracy of iterating V X E to zero. 

The validity of the PRECHG CODE has been checked numerous times by comparison 
of the results of calculations with simple geometries that can be solved in 
closed form and by comparison to the results of other analytical computer 
codes ouch as the Arbitrary Body of Revolution Code (ABORC) which themselves 
have been previously verified. 
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